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Using convex optimization, we propose entanglement-assisted quantum error correction procedures that are 
optimized for given noise channels. We demonstrate through numerical examples that such an optimized error 
correction method achieves higher channel fidelities than existing methods. This improved performance, which 
leads to perfect error correction for a larger class of error channels, is interpreted in at least some cases by 
quantum teleportation, but for general channels this interpretation does not hold. 
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I. INTRODUCTION 

Noise is an important obstacle for the scale-up of quan- 
tum information devices — both decoherence due to interac- 
tions with an external environment, and the internal control 
errors inevitable for any information-processing. The theory 
of quantum error correction was developed by analogy to clas- 
sical error-correcting codes to overcome this obstacle [1—4]. 
Quantum error-correcting codes (QECCs) store quantum in- 
formation redundantly in an entangled state of multiple quan- 
tum systems (usually qubits), in such a way that errors can be 
detected and corrected without directly measuring (and hence 
disturbing) the quantum information to be protected. The fact 
that this can be done at all is quite remarkable, and has opened 
up the prospect that quantum information processing may be- 
come a viable technology in the foreseeable future. 

Broadly speaking, quantum error correction is useful in two 
main contexts. The first is quantum computation, where errors 
that occur during a calculation must be constantly found and 
corrected to keep them from accumulating. The QECCs used 
for this purpose must fit into an overall fault-tolerant scheme 
[5-8], which puts important constraints on the properties of 
the codes used. 

The other main context is quantum communication. Here it 
is assumed that the steps of encoding and decoding are largely 
error-free, and all errors occurs during the transmission of 
quantum information through a noisy quantum channel. This 
is the context we will consider in this paper. It should be 
noted, though, that this division is not really sharp: quantum 
computers may well need to transmit information internally, 
and quantum communication may be part of a distributed in- 
formation processing scheme. 

In the context of quantum communications, one can ex- 
plore the effect of different shared resources on the ability 
to detect and correct errors. Recently, a great deal of work 
has been done on entanglement-assisted quantum error cor- 
rection (EAQEC), and entanglement-assisted quantum error- 
correcting codes (EAQECCs). In these schemes, the sender 
(usually called Alice) and receiver (usually called Bob) share 
some number of pure maximally-entangled states, or ebits. 



Alice can use her halves of these ebits in the encoding pro- 
cedure, and Bob can use his halves in error correction and 
decoding. The idea of using entanglement to improve error 
correction was proposed early in the history of quantum in- 
formation theory [9], but it took a long time for the first steps 
towards designing practical codes [ ]. This work has led to 
a generalized theory of quantum error correcting codes, the 
entanglement-assisted stabilizer formalism [1 1-17]. 

The large majority of work on QECCs has concerned sta- 
bilizer codes [ ], which can be derived from classical linear 
codes. In this algebraic approach, one generally tries to de- 
sign codes that are able to detect and correct an abstract set of 
errors acting on single qubits, often based on the Pauli opera- 
tors. QECCs have the property that any linear combination of 
correctable errors is also a correctable error. Since the Pauli 
operators form a basis for all single-qubit operators, the abil- 
ity to correct Pauli-based errors implies the ability to correct 
more realistic models of noise, provided that errors are not too 
highly correlated between qubits. The EAQECCs that have 
mostly been derived so far extend this stabilizer formalism to 
include a broader class of codes that utilize shared entangle- 
ment. 

The standard approach to QEC finds encoding and recovery 
procedures which ensure perfect recovery of quantum states 
passing through sufficiently weak noisy channels. The dis- 
advantage of this approach is that it may fail when the noise 
strength, or error probability, is too high. Another approach 
to QEC is to tailor QECCs to a particular model of realis- 
tic noise, derived from a detailed model or from experimental 
measurements [ 1 8-20], e.g., via quantum process tomography 
[21-27]. The tailored approach to QEC attempts to find en- 
coding and recovery operations which are optimized relative 
to the particular experimentally measured noise or assumed 
noise model, in the sense that they maximize the fidelity (or 
minimize the distance) of the encoded output state relative to 
the input state [28-39]. Both worst case and average case per- 
formance has been considered in this setting, which usually 
involves numerical optimization. The advantage of this ap- 
proach is that it tends to be more robust to variations in the 
noise strength than the standard approach [35]. Here we fur- 
ther develop the optimized QEC approach by considering, for 
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the first time, entanglement-assisted quantum error correction 
as an optimization problem. 

To formulate the problem for optimization, we break down 
quantum error correction into three stages. First, the quantum 
state to be protected is encoded. This procedure appends some 
number of ancilla qubits in a fixed initial state, and then ap- 
plies an encoding unitary to the information and ancilla qubits 
together. These qubits then pass through a noisy channel (as- 
sumed to be known). At the receiver's side, more ancillas can 
be appended, and then a decoding unitary is applied. All the 
ancillas (encoding or recovery) are discarded, and the state of 
the information qubits is the output. 

In standard quantum error correction, the encoding and re- 
covery ancillas always start in a standard state (usually |0)), 
and are not entangled with each other. Previous work has 
shown that regular recovery ancillas are redundant in increas- 
ing the error correction fidelity [35]. That is, given the same 
number of information qubits and encoding ancillas, the pres- 
ence or absence of recovery ancillas makes no difference to 
the channel fidelity. This matches the properties of stabilizer 
codes, where decoding can be done unitarily, and there is no 
benefit in including recovery ancillas. 

To formulate the problem for entanglement-assisted codes, 
we change the initial state of the ancillas, so that noiseless en- 
tanglement is shared between the encoding and recovery an- 
cilla qubits. We show that this additional information can en- 
hance the performance of the recovery ancillas, and increase 
the error correction fidelity in most cases. This initial shared 
entanglement can in principle improve performance in at least 
two different ways. It can allow the code to transfer quantum 
information about the initial state to the recovery block, im- 
proving the channel fidelity. This fidelity increase can in fact 
lead to perfect error correction for an important class of er- 
ror channels. The error correction procedure for this class of 
channels can be interpreted as teleportation, with only classi- 
cal information passing through the noisy channel. 

However, even in cases where the error correction protocol 
is not a form of teleportation, we can still get an improvement 
in performance. One interpretation of this is that entangled 
ancillas allow the recovery procedure to extract more infor- 
mation about the errors. This is analogous to superdense cod- 
ing, where the use of entanglement boosts the rate of classical 
communication. The fidelity increase here shows the impor- 
tance of relevant information stored in the ancillas in quantum 
error correction. The entangled recovery ancillas are funda- 
mentally different than the regular recovery ancillas consid- 
ered in previous optimization problems, that gave no benefit 
[35]. 

The optimization procedure here, similar to previous work 
in this area [28-39], assumes that we know the noise chan- 
nel. We assume that a channel identification procedure, such 
as quantum process tomography [ 1 8-27], has been performed 
prior to the error correction. It is also possible to optimize 
performance in cases where the exact channel model is not 
known, though we do not consider that here. 

The structure of the paper is as follows. In Section II we 
formulate the optimization problem in terms of a distance be- 
tween states that needs to be minimized. The problem is bi- 



convex in the encoding and recovery operations, and in Sec- 
tions III and IV we consider how to optimize the encoding for 
a given recovery, and recovery for a given encoding, respec- 
tively. In Section V we consider the case of random unitary 
channels, and analyze a simple example for which EAQEC 
has perfect fidelity, via teleportation. We present numerical 
examples in Section VI, where we contrast different correc- 
tion and optimization scenarios. These examples illustrate 
the improved performance of optimized EAQECCs in a set- 
ting where perfect fidelity cannot be achieved. We conclude 
in Section VII. Certain technical details are presented in Ap- 
pendix A. 

II. PROBLEM FORMULATION 

A quantum dynamical process is a map from an initial state 
p to a final state £{p). In general the map £ does not repre- 
sent unitary evolution. However, if we assume a larger closed 
quantum system that includes the environment (or bath) B of 
the system, the total evolution becomes unitary. If the initial 
system-bath state is psb, a quantum channel can be written as 

£(p)=Ttb[U S bPsbUI b }, (1) 

where p — Tr# [pg g] is the initial system density matrix, Us b 
is a unitary operator acting on the joint system-bath Hilbert 
space % — Hs^T-Lb, Xb — ^s[psb] is the initial state of the 
environment, and Tr^ (Trg) denotes the partial trace over the 
environment (system). In general a quantum dynamical pro- 
cess described by Eq. (1) is a Hermitian map (i.e., it maps Her- 
mitian operators to Hermitian operators) [40]. When the ini- 
tial system-bath state is purely classicaly correlated (formally, 
has vanishing quantum discord [41]), f becomes a completely 
positive (CP), trace-preserving map [ ]. This condition is 
also necessary [ ]. In this case the evolution can be rep- 
resented using the Kraus operator sum representation (OSR) 
[43^4-5]: 

£{p)=Y,Ki P Kl (2) 

i 

where the operators Ki, known as Kraus operators, satisfy the 
normalization condition J^i — I (identity). Such maps 

are powerful tools which can represent different types of noise 
processes that could be acting on the system. 

We shall assume that the error correction procedures of en- 
coding and recovery, as well as the noisy channel, can all be 
described using Kraus operators. The entire procedure on the 
system can then be represented as 

C £ 1Z 

Po Pc Pe pr, (3) 

where C, £, and 1Z are the encoding, error and recovery op- 
erations respectively. The density matrices po through pn are 
the states of the data qubits and the ancillas, i.e., represent the 
state of the entire system at different points in time. The ini- 
tial system state p : Hs where Hs is the total system 
Hilbert space, is taken to be a product state between the data 
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qubits, represented by pdat, and the ancillas, represented by 
Pane, i-e., po = Pdat ® Pane- Thus % s = H da t ® "H a nc, where 
%dat and ^anc are, respectively, the Hilbert space of the data 
qubits and of the ancillas. The ancilla Hilbert space further 
decomposes into an encoding subspace H cnc and a recovery 
subspace W. rcc : "H an c — 'H cnc ®'Hrec- We write d = dim('Hs) 
and otherwise denote the respective dimensions of the various 
subspaces by: 

dx — dim('Hx), X £ {dat, anc, enc, rec}. (4) 
Note that 

d ^dat^anc; d anc — d cnc di Cc . (5) 

For most of our discussion we shall assume that all encod- 
ing ancillas are pairwise maximally entangled with recovery 
ancillas prior to the encoding operation, whence d cnc = d ICC , 
and p anc — (\B) (£>|)® <ionc , where \B) is a maximally entan- 
gled state between encoding and recovery ancilla pairs. Nev- 
ertheless, for generality we shall leave d cnc and d rcc unless it 
is necessary to be more specific. Later, in our numerical ex- 
amples, we shall also consider regular (unentangled) encoding 
and recovery ancillas. 

Using the OSR for each of these operations, the overall evo- 
lution is 

p R = TZ o £ o C(po) 

= Y,( R r E eCc)po(RrE e C c )\ (6) 

c,e,r 

where {Cc}™^, {E e } r e %, and {R r }T" x are the Kraus oper- 
ators of the encoding, error, and recovery channels, respec- 
tively, all satisfying the normalization condition: 

C\C C = E l E - = R r R r = L ( 7 ) 

c e r 

All the Kraus operators are here represented by square d x 
d matrices, i.e., are operators on Hs- The number of Kraus 
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FIG. 1: (Color online) Scheme of entanglement-assisted error cor- 
rection. The data qubits are represented by pdat • The encoding and 
recovery ancillas are entangled before the encoding procedure C be- 
gins, and are represented by p a nc- The joint initial system state is 
Po = Pdat <8> Pane- The recovery ancillas pass through a noiseless 
channel to the recovery operation TZ. The information qubits and en- 
coding ancillas pass through a noisy channel £ . The output of the 
recovery procedure is pa- The intermediate states pc (output of the 
encoding) and pe (output of the noise channel) are also indicated. 



operators for each map depends on the implementation and 
the basis representation, but the number txie of error operators 
is bounded above by d 2 [44]. 

We shall assume that the encoding operation C is unitary, 
represented by a single unitary (Kraus) operator C. The en- 
coded d-dimensional state is therefore 

PC =^(p da t®p anc )C t . (8) 

Since, as illustrated in Figure 1 , the encoding does not act on 
the second half of the entangled pairs (the recovery ancillas), 
the encoding unitary C should be represented as C ® ircc, 
where C is a dc x dc matrix, and where dc> = rfdat^enc = 
ddat Vd&nc- We shall write Ix to denote the identity operator 
acting on a subspace X, and reserve / for the identity operator 
on the entire system space Hs- 

Our goal is to implement a close approximation to a certain 
desired unitary L on the entire system space. More specifi- 
cally, the purpose of optimization is to design the encoding 
C and recovery TZ for a given error channel £ such that the 
overall map (6) po H> pb. is as close as possible to some de- 
sired unitary L : T~Ls > T~Ls- This goal is different from 
that pursued in earlier optimization papers (e.g., [33-35]), 
where the goal was to obtain a good approximation to a uni- 
tary Ldat : "Hdat | -> %dat acting only on the data qubits. Our 
present goal is general enough to encapsulate state preserva- 
tion of the data qubits, e.g., by implementing, with as high a 
fidelity as possible, a unitary which swaps the data qubits into 
the recovery qubits at the end of the process, while leaving 
the data and encoding qubits in a new pure state which is ten- 
sored with the recovery qubits. Indeed, this is essentially what 
happens in our teleportation example (Section V). 

Channel fidelity or distance are typical measures of perfor- 
mance between two quantum channels [28, 44, 46, 47]. The 
channel fidelity between the error correction operation 1Z£C 
and the desired unitary operation L is: 

f = ^J2\ TrL ^ R r E e C \ 2 > (9) 

We wish to maximize this fidelity, which satisfies < / < 1. 
From [ , Thm.8.2], / = 1 if and only if there are constants 
p re such that 

R r E e C = p re L, ^| Mre | 2 = l. (10) 

As shown in [35], maximizing the fidelity / is equivalent to 
minimizing the following distance, which is the approach we 
focus on henceforth: 

S = J2\\ R r E eC- PreL\\ 2 , (11) 
r,e 

where ||X|| = (TyX^X) 1 / 2 is the Frobenius norm. The mini- 
mization must now include the parameters {p re } as well. Us- 
ing straightforward matrix manipulations, we can rewrite this 
distance as 

S(R, C, A) = \\RE(I rnE «> O - A <g> L\\ 2 (12) 
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where A = [/x re ] is an run x uie rectangular matrix, E is the 
d x niE<i rectangular error matrix built using the error Kraus 
operators as horizontally arranged blocks 



E — [E\...E mE \, 



(13) 



and R is the mud x d rectangular matrix obtained by vertically 
stacking the recovery Kraus operators R r : 



R 



Rn 



(14) 



Therefore, we have || A|| 2 = Tr A^A = e \P"re\ — 1> an d 

R)R = J2 r RtRr = I- 

Our optimization problem can now be summarized as fol- 
lows: 

minimize 5{R 1 C, A) subject to the constraints 

l|A|| 2 = l, 
(15) 



rIr = I, C = C'®I Iec , C'^C = I t 



dat.: 



where idat,rec denotes the identity operator on Hd&t <8> ^rec- 
There are three matrices of parameters (R, C, A) in this ex- 
pression that are to be optimized. As stated, this optimization 
is not convex. Therefore, as first suggested by Reimpell and 
Werner [28], to identify the optimal values of these param- 
eters, we solve this optimization problem iteratively. Start- 
ing with an arbitrary initial encoding operation, we first find 
the optimized recovery operation, which by itself is a con- 
vex problem; we then fix the recovery operation obtained in 
the previous step, and find the optimized encoding operation, 
which is also a convex problem; and so on, alternately opti- 
mizing the recovery and encoding operations until the proce- 
dure converges. A is an intermediary matrix of parameters 
that must also be recalculated at each step. This iteration con- 
tinues until the distance (11) stops decreasing. The details of 
the procedure at each step are provided in the next two sec- 
tions. 

The design of the optimized recovery operation at each step 
is not affected by the extra constraints added by introducing 
entanglement in the error correction procedure. This is so be- 
cause optimization of the recovery operation does not depend 
on the details of the encoding operation. However, the ex- 
tra constraints play an important role in identifying the opti- 
mized encoding operation, which must be of the desired form 

C — — C I^ec ' 

The encoding ancillas must be initialized in a maximally 
entangled state — for example, an EPR pair. However, the 
distance measure defined in (11) is independent of the initial 
state of the ancilla qubits. To overcome this problem we add 
an extra step, an entangling operation U : Hs H> Hs, to 
the procedure. The entangling operation does not act on the 
data qubits, so it can be written as U = /dat <£> U anc , where 
C^anc : Ham >-> H a nc- By including this entangling operation 
in the evolution, the ancillas can be assumed to have the initial 
state |0). 



III. OPTIMIZED ENCODING OPERATION 

The derivation of the optimized encoding operator is simi- 
lar to what was done previously in [ 55], but requires includ- 
ing the entangling operator U, and a careful consideration of 
the subspaces involved. Including the entangling operator, the 
distance measure (11) can be rewritten as 

6{R,C,A) = \\RE(I mE ® CU) - A ® L\\ 2 

= \\RrE e CU - H re Lf 

r,e 

= Tr[I + C^C-UL^Y,^reRrE e C-h.c.}, 

r,e 

(16) 

where in the second line we used the relation between 
Eqs. (11) and (12), and in the third line we used the defini- 
tion of the Frobenius norm and the normalization conditions 
(7) and (10). 

As a first step we'll need to find the unconstrained mini- 
mum of S(R, C, A) with respect to C, i.e., we need to find 
the solution to 85 / dC = without introducing the condition 
C^C = I. To this end we note the following matrix differ- 
entiation identities, valid for any pair of matrices A and Z 
with compatible dimensions, with A being independent of the 
elements of Z: 



d_ 
~dZ 
d_ 
~dZ 



TvlAZ] 



Tr[ZZ f ] 



.4' 



Z*, 



d_ 
~dZ 
d 



dZ 



Tv[A{Z (g> I)} 



dZ 
(Tr 2 A)\ 



Tr[AZ^} = 0. (17) 
Tr[AZ^Z] = Z*A\18) 
(19) 



where in the last identity the partial trace is over the subspace 
acted on by the identity operator. We prove all these identities 
in Appendix A. 

Therefore, assuming that the recovery operators R r are 
fixed, we have 



95 
dC 



C* 



,RrE e ) 



(20) 



so that the solution to the unconstrained optimization problem 

— — is 
dc — u 1!> 



C = J2Vre{R r E e ^LU 



(21) 



We wish to impose the additional structure C — C'cg/roc, i.e., 
C = ~r— Tr rec C, where the partial trace is over the recovery 
ancillas. Thus the unconstrained solution over "Hdat <8> %enc 
becomes 



C = argmin(5(i?,C'<8>/ r oc,A) 



-Tr r 



J2Vre(R r E e yL(I dat ®Ul) 



.(22) 
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The reader can check that this result can also be derived by 
introducing the decomposition C — C ® / rec directly into 
Eq. (16) and using identity (19), i.e., 

88 
dC' 



d iec C* - (Tr nc [ULlY,6e R r E e]) t - (23) 



This is important for the next step. 

The constrained problem (15) requires that C'^C = 
7dat,rec- To solve this problem we form the Lagrangian 

C = S(R, C, A) + Tr[P(C"tC" - 7 dat>rec )], (24) 

where P is a Hermitian Lagrange multiplier matrix (be- 
cause the constraint C'^C = /<iat,rec is Hermitian). Set- 
ting dc"C = we find, using Eqs. (18), (22), and (23): 
C* - C* + C"*P' = 0, i.e., the solution to the constrained 
problem is 

C = C'{I dat ^ c + P)-\ (25) 

To eliminate P we note that the constraint C'C = /dat.rec 
now implies C*C' = (ldat,rec + P) 2 , from which we have 
(/dat,rec + P)' 1 = (C'W)- 1 ' 2 - Thus, finally 



(26) 



where C is given in Eq. (22). This encoding operation is 
optimized for the given recovery operation R. 



IV. OPTIMIZED RECOVERY OPERATION 

Unlike the case of the optimized encoding operation, we do 
not impose any specific tensor product structure on the recov- 
ery, so that derivation of the optimized recovery operator is 
essentially the same as what was done previously in [35], ex- 
cept that again we must include the entangling operator. Here 
we briefly outline the steps required to optimize the recovery 
operator. Using the constraints in (15) and fixing the encoding 
operation C, we express the distance (1 1) as 



S = \\RE(I mE ® CU) - A ® L\\ 
= 2d-2ReTri? J B(A t «)CC/L t ). 



(27) 



Minimizing (27) with respect to R is the same as maximiz- 
ing the last term. As shown in [ 15, Appendix A] this maxi- 
mization yields 



max Re Tr REfA^ <g) CUtf ) 

&R=I 



= Try E(T ® J)Et 
where the matrix T is defined as 

r = A f a. 



(28) 



(29) 



The constraint || A|| =1 from (15) is equivalent to Tr T = 
1, and r > by definition. Therefore, the optimization prob- 
lem for T is: 



maximize Try' E(T ® I) Eft subject to 
T > 0, TiT = 1. (30) 



The optimized T can be obtained by solving an equivalent 
semidefinite programming (SDP) problem, or by performing 
a constrained least squares optimization [35]. The matrix A 
can be obtained from T by the definition (29) (see Eq. (21) 
of Ref. [ ] for details). Once the optimized A is known, the 
optimized recovery matrix can be found from 



R = [vi ... v d ][ui ... u d ] 



(31) 



where {vi,Ui}f =1 are, respectively, the right and left singu- 
lar vectors in the singular value decomposition of the matrix 
E(A* <g> CU U), with the singular values in descending order. 

The algorithm below summarizes the preceding method for 
encoding and recovery optimization in the presence of entan- 
glement between the encoding and recovery blocks. 

Initialize C 
Repeat 

a) Optimal recovery 

maximize (30) to find T 
solve for A via (29) 
solve for R via (31) 

b) Optimal encoding 

solve for C via (26) 
Until distance (11) stops decreasing 

Since in each iteration of this algorithm the distance 5 only 
decreases, the converged solution to this optimization is guar- 
anteed to be at least a locally optimal solution to (15). We 
apply this procedure to a selection of different channels below 
to assess its performance in practice. 



V. RANDOM UNITARY CHANNELS 

An error channel £ is called a random unitary channel if 
we can decompose it into the probabilistic application of one 
of a finite set of unitary operations: 



£(f)=Yjp i V i pV! ) 



(32) 



where the Vi are unitary operators and p^,i = 1, . . . , n is a 
probability distribution. Random unitary channels describe 
noise processes that can be corrected using classical infor- 
mation extracted from the environment [ ]. It is clear that 
any random unitary channel should be unital, meaning that 
£ (I) = I. However, the inverse relation holds only for chan- 
nels on qubits, and is not true for higher dimensions [49, 50]. 

There are two operations in ( 1 ) that are not random unitary. 
These operations are the introduction of the ancillas and the 
partial trace over the environment. Therefore, an error chan- 
nel need not be random unitary, of form (32), in general. The 
necessary and sufficient conditions for a channel to be a ran- 
dom unitary are discussed in [51]. In [52], an upper bound 
was found on the number of unitaries needed in (32). Here 
we show that if the number of unitary operators required in 
(32) is at most two, our optimized error correction method can 
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perfectly correct the error using only one EPR pair as ancilla 
qubits. 

In fact, for this case, the optimized encoding and recovery 
operations are known, and are equivalent to teleportation. To 
see this, consider a channel that can be decomposed into two 
unitaries, V\ and V 2 , as follows: 



£{p) = {l-p)V lP Vl + P V 2 pV, 



rt 



(33) 



Suppose Alice wants to send Bob one qubit of data through 
this noisy channel. One pair of entangled qubits enables Alice 
and Bob to make this communication error-free. Assume this 
pair of qubits was prepared by an entanglement source and 
sent to Alice and Bob, each taking one of the qubits. Alice 
can apply an encoding operation on her qubits — the data qubit 
and her half of the ebit — before sending them to Bob through 
the noisy channel. Bob applies the recovery operation to all 
three qubits once he receives the message. 

The optimized encoding operation that Alice uses in this 
scenario is done by first applying V± to her qubits (recall that 
we assume that the channel is known), to make the error chan- 
nel equivalent to 

£{p) = (1 - p)I P I + pV^VlpVxVl (34) 
Alice continues the encoding operation by applying 

C = ^(\v 1 v 1 )(B 1 \ + \v 1 v 2 )(B 2 \ 

+ \v 2 v 1 )(B 3 \ + \v 2 v 2 )(B 4 \), (35) 

where v\ and v 2 are the eigenvectors of the unitary V 2 Vl, 
and B\ to B4 are an orthonormal basis of maximally entan- 
gled states — for example, the Bell states (|00) ± 1 1 1) ) /\/2 
and ( 1 10) ± 1 1) ) /V2. This encoding enables us to send two 
bits of classical information safely through the noisy channel. 
These two bits of classical information, together with the sec- 
ond half of the ebit, enable us to recover the original state. 

To see this in more detail suppose that the initial state of 
the data qubit is \(f>i) — a |0) + b The state of the entire 



l?o> 



|9i> [H_ 
W2) 



e 



ns 



x - z - 



FIG. 2: Quantum error correction using teleportation. The infor- 
mation qubit \qo) is recovered with unity fidelity in spite of passing 
through the noise channel £ . The first two operations (Hadamard 
gate on |qi) and controlled-NOT between \qi) and | Q2 } ) represent 
the entangling operation. This is followed by the encoding operation 
C described by Eq. (35). The last two operations represent condi- 
tional X and Z gates applied to \ q2), dependent on the measurement 
outcomes of |go) and |gi). See the text for additional details. 



system after applying the entangling unitary is 

hfc) = ^|K|000) + |011>) + 6(|100) + |111»]. (36) 

After applying the encoding (35) and error map on the first 
two qubits, the final state of the entire system is 



1 



(|«i«i)(o|0) + 6|l)) 



2V2 

+ \v 1 v 2 )(a\0)-b\l)) 
+ \v 2 v 1 )(a\l)+b\0)) 
+ \v 2 v 2 )(a\l)-b\0))). 



(37) 



This is the state of the system once Bob receives it. 

At this point, as shown in Figure 2, he can recover the orig- 
inal state \4>i) by measuring the state of the first two qubits. 
Based on the result of the measurement he applies the appro- 
priate Pauli operator to the last qubit, his half of the ebit, to 
recover the original state. For example, if he measures the 
first two qubits to be \v 2 v\), he should apply a x to the last 
qubit to recover the original state. 



VI. EXAMPLES 

A. The Bit Flip Channel 

As an example of a random unitary channel with two uni- 
tary operators, we consider the bit flip error channel, or bi- 
nary symmetric channel, in which bit flip errors occur inde- 
pendently with probability p on each qubit. This channel can 
be represented as 



£{p) = (1 - p)lpl + pa x pa x , 



(38) 



where henceforth a x ,a y , a z denote the standard 2x2 Pauli 
matrices. As shown in Figure 3, this error can be corrected us- 
ing an ebit as the error-correcting resource. The ebit again is 
shared between the encoding and recovery parts. The encod- 
ing acts on the initial data qubit and the first entangled qubit 
(the encoding qubit). Considering that the error is already in 
the form (34), we do not need to apply the initial unitary oper- 
ator here. Therefore, the optimized encoding operator based 
on (35) is 



C = 



(Si I 



-) (B 2 



(B 3 \ 



I—) (Ba 



(39) 

where |+) = (|0) + |1»/V2, R = (|0)-|1»/V2, and Si to 
_B 4 form an orthonormal basis of maximally entangled states. 
The first entangled qubit together with the initial data qubit 
are corrupted by the bit flip error while the other entangled 
qubit (the recovery qubit) stays intact. The recovery operation 
acts on all the qubits, and can reproduce the initial state by 
measuring the first two qubits as discussed above, followed 
by a SWAP between the recovery qubit and the data qubit, 
which perfectly restores the state of the data qubit. 

While the optimization procedure we employ is not guar- 
anteed to find the global maximum of the channel fidelity, we 
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see that in this case it does indeed find the optimal solution. 
Convex optimization has rediscovered the protocol of quan- 
tum teleportation. 



B. The Bit Flip/Phase Flip Channel 

In this example both bit flip and phase flip errors occur with 
equal probabilities: 

P P 

£{p) = (1 -p)lpl + -<r x pa x + -<J z p<r z . (40) 

The error occurs independently on different qubits. This er- 
ror channel clearly cannot be represented in the form of (33). 
Hence our teleportation argument does not apply here. There 
is no choice of basis for which classical information can be 
sent through this channel error-free. 

While entanglement cannot enable perfect error correc- 
tion in this case, it can still increase the useful information 
about the errors available to the recovery block, and there- 
fore increase the fidelity of the error correction. Figure 4 
presents the fidelity of error correction for this channel for 
different values of p. The only difference between the up- 
per and middle curve is the presence of entanglement. With- 
out entanglement, the extra qubit used in the recovery acts 
as a regular recovery ancilla, and therefore does not increase 
the fidelity. The minimum fidelity in both cases occurs at 




0.2 0.4 0.6 0.8 1 

P 



FIG. 3: (Color online) Optimized fidelity for the bit flip channel. Re- 
sults are shown for the case of one data qubit. The blue (bottom) 
curve is the fidelity without any recovery. The red (middle) curve is 
the optimized fidelity obtained for encoding using a single encoding 
qubit but for recovery without entanglement, i.e., with a single recov- 
ery qubit in a product state with the encoding qubit. This represents 
the standard scenario of (channel-optimized) QEC. The green (top) 
curve is the fidelity obtained in the presence of a single ebit, i.e., the 
EAQEC scheme shown in Fig. 1 . The fidelity in this case is unity for 
all values of the bit flip probability p, showing that our optimization 
recovers the perfect encoding and recovery protocol of teleportation 
discussed in the text. 



p = 2/3, where the channel is symmetric in /, a x and a z \ 

£{p) = (l/3)(lpl + <7 x p<7 x + a z pa z ). 



C. The Depolarizing Channel 

In this example we consider a particularly simple error 
model that is widely used to describe noisy quantum systems: 
the depolarizing channel. This important channel is often con- 
sidered to be the quantum equivalent of the binary symmetric 
channel in classical error correction. It occurs when the state 
of the system can be completely mixed by the action of the 
channel. The noise is unbiased in the sense that it generates 
bit flip errors, phase flip errors, or both with equal probability; 
this is represented by the map 

P 

£{p) = (1 - p)p + - (<J x pcr x + Oypoy + a z pa z ) . (41) 

The depolarizing channel shrinks the radius of the Bloch 
sphere by a factor of 1 — p, while preserving its shape [44]. 

Interestingly, for this error channel the ebit does not in- 
crease the fidelity for p < 3/4. The optimized fidelity of 
the channel with and without using an ebit is presented in Fig- 
ure 5. A break point occurs at p = 3/4, when the coefficients 
in (41) become equal. The fidelities of both cases are the same 
for p < 3/4, but the ebit can slightly increase the fidelity for 
p>3/4. 

An observation that may explain this difference is that for 
p < 3/4 we can write £{p) = (1 - 4p/3)p + (4p/3)T(p), 
where T is the twirling operation: 

T(p) = (p + cr x pa x + o y pay + a z pa z ) = -I, (42) 



1 r 




P 



FIG. 4: (Color online) Optimized fidelity for the bit-phase flip chan- 
nel. As in Fig. 3 there is a single data qubit going through the noise 
channel unprotected (blue, bottom curve), optimally protected by a 
single encoding and a single recovery qubit in a product state (red, 
middle curve), and optimally protected by an ebit (green, top curve). 
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for any valid density matrix p. In other words, with probability 
4p/3 the channel replaces the qubit with a maximally mixed 
qubit. For p > 3/4 this interpretation is no longer possible, 
and it is in this regime that shared entanglement seems to give 
improvement. 

VII. CONCLUSION 

Previous work on optimized quantum error correction [28- 
39]. considered the "standard", entanglement-unassisted set- 
ting. In this work we addressed, for the first time, the problem 
of optimized entanglement-assisted quantum error correction. 
Our methodology is a generalization of the optimization ap- 
proach developed in Ref. [ ], which involves a bi-convex op- 
timization problem, iterating between encoding and recovery. 
Our main technical innovation appears in the encoding step, 
where we introduced a method to account for the constraint of 
initial entanglement between the encoding and recovery ancil- 
las. 

Our method returns the optimized channel fidelity for a 
given noise channel in the presence of shared entanglement 
between the encoding and recovery block. We showed that 
entanglement can substantially increase the error correction 
fidelity. For specific error channels this fidelity increase can 
lead to perfect error correction. In these cases, the error cor- 
rection procedure can be interpreted as quantum teleportation. 
For cases in which perfect error correction is not possible, this 
teleportation interpretation does not apply; however, our opti- 
mization method returns the optimized fidelity for these chan- 
nels, in most cases showing improvement with the resource of 
shared entanglement. 

For one channel we found no improvement: the depolariz- 



0.9 r 




P 



FIG. 5: (Color online) Optimized fidelity for the depolarizing chan- 
nel. As in Figs. 3 and 4, there is a single data qubit going through the 
noise channel unprotected (blue, bottom curve), optimally protected 
by a single encoding and a single recovery qubit in a product state 
(blue, middle curve, overlapping red for p < 3/4), and optimally 
protected by an ebit (red, top curve). 



ing channel. For the small codes that we considered in this pa- 
per, optimized standard quantum error correction performed 
just as well as optimized entanglement-assisted quantum er- 
ror correction on this channel. It is possible that this equiv- 
alence between entanglement-assisted and standard quantum 
error correction will not hold for larger codewords on the de- 
polarizing channel. Larger codewords would use either more 
ancillas or more ebits in the encoding process. The possibility 
that larger codes might show a greater difference is suggested 
by the fact that a three-qubit codeword with two ebits exists 
that can correct an arbitrary single-qubit error [ ], while the 
smallest standard code that corrects an arbitrary error is of size 
five [53]. However, the performance of optimized quantum 
error correction for larger codewords remains a subject for fu- 
ture research. Larger codewords make for a more numerically 
intensive optimization procedure. 

Considering that regular recovery ancillas are redundant in 
the optimized error correction procedure, the results of this 
paper show that information about errors in the encoded data, 
transferred to the recovery block through entanglement with 
the encoded state, can increase the fidelity for many channels. 
We expect that results for a wider variety of error channels 
and larger codewords will continue to bear this out. 
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Appendix A: Trace differentiation identities 

We prove Eqs. (17) and (18). These are one-line proofs 
which follow directly from the definition of differentiation by 
a matrix [54, Section 10.6]: the derivative -J^f(Z) of a dif- 
ferentiate scalar-valued function f(Z) of a matrix argument 
Z = [zij] E c« ix ™ i s the m x n matrix whose entry is 

dz Tl[AZ] = = M = A *> < A1 > 

%1 ij 

±Tr[A2fi] = [^1^41=0. (A2) 

13 ij 

^Tr[AZ^Z] = [■^ y £z iS (Atf) ji ] = [{Atf) ji ] 

13 ij 

= Z*A\ (A4) 

where for the second identity we used the fact that for a com- 
plex variable z = x + iy (x, y E W) it holds that dz* jdz — 
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d(x — iy)/dx — id(x — iy)/dy = 0. 
Next we prove Eq. (19). Let 

A = ^a m n P q\m){n\ ® \p){q\, (A5) 
Z = Z®I = ^z mln ,\m'){n'\®\p')(j>'\. (A6) 

Then 

Tr AZ = Tr ^ a mnpp > z nn , \m){n'\® \p) {p' \ 

— ^ y Q'mnppZnmi (A7) 
mnp 



so that 

d d 
—Tr[A(Z (g> I)} = [- ^2 a mnp pz nm } 

aA " Zmn mnp 

= [Y, a nm P p] = (Tr 2 A) { , (A8) 
p 

where the partial trace is over the subspace acted on by the 
identity operator in Z ® 7. 
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